Efficient readout of flux qubits at degeneracy 
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We present a superconducting circuit consisting of a flux qubit and a single-charge transistor 
serving as a detector. As flux and charge are conjugate, the transistor can detect states of the qubit 
close to the flux degeneracy point, when the eigenstates are quantum superpositions of fluxes. The 
coupling has a flip-flop symmetry conserving the total number of excitations, and so the measure- 
ment outcome results in the absence or presence of an incoherent tunneling cycle. We evaluate 
the performance of a practical device and show that it is an attractive tool for measuring at the 
degeneracy point. 
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Small Josephson junction circuits are a prototype sys- 
tem for studying generic quantum effects in solid-state 
physicsi These devices are a candidate for scalable 
quantum computing, which requires full and precise con- 
trol of quantum states in real time. The strong devel- 
opment in this field has lead to a number of substan- 
tial achievements such as superpositions^ and real-time 
control^ of collective states as well as the control of charge 
states^ and the demonstration of a controlled-not gate^. 
A highly attractive realization of a superconducting qubit 
is based on magnetic flux states^. 

A crucial and subtle issue in this effort is the proper 
understanding and design of the read-out device. Quan- 
tum measurements are intrinsically invasive^. Ideally, the 
back-action exclusively consists of instantaneous projec- 
tion on eigenstates of the measured observable. In gen- 
eral, and specifically for quantum systems such as the 
Josephson circuits, the back-action consists of two parts: 
The dcphasing process which projects onto the measure- 
ment basis, and relaxation which induces transitions be- 
tween these states leading to errors in the observed vari- 
able. The former is inevitable and is an integral part of 
quantum control protocols, whereas the latter has to be 
avoided^. This aspect makes quantum nondcmolition 
(QND) schemes highly attractive. Thus, for fundamen- 
tal as well as for practical reasons we are faced with a 
problem of understanding and ultimately engineering the 
properties of the read-out device on the quantum level. 

Next to the backaction of the detector, any quantum 
bit realized in a solid-state environment suffers from de- 
coherence induced by the other degrees of freedom of 
the solid-state environment and of the fluctuations of the 
control fields — flux noise in the case of flux qubits. Slow 
fluctuations from the ubiquitous 1// noise in the solid- 
state environment are particularly destructive, here re- 
sulting mostly in critical current fluctuations. Thus, it is 
very attractive to operate devices at the charge degener- 
acy point where they are to first order immune to these 
fluctuations&iLiS. It is possible to achieve a very stable 
intrinsic bias at that point by flux trappingifi. Fixing 



the bias to that point can lead to restrictions on the con- 
trol parameters. Although operation by AC-pulses is still 
possible, it is not evident how to read out the state at 
the end. In a flux-based device operated at flux degener- 
acy any Hamiltonian-dominatecU*S measurement cannot 
distinguish the energy eigenstates of the system. The 
natural solution is the measurement of a different vari- 
able such as phased, which in the particular case of Ref. 
|3l9l | still requires to move out of the optimum point. 

In this letter, we are going to show how, using a 
Cooper-pair transistor, the state of the qubit can be read 
out with high fidelity and low unwanted backaction with- 
out leaving the degeneracy point. The physical coupling 
mechanism is based on geometric phases. By appropriate 
biasing of the transistor and by using a readout chan- 
nel which filters out unwanted processes, the physical 
processes are flip-flop transitions between qubit and de- 
tector, which lead to a QND-like and quantum limited 
backaction. We will discuss values of experimentally ac- 
cessible parameters. 

At the degeneracy point, the energy eigenstates are 
coherent superpositions of flux states, which cannot be 
discriminated by flux detection. On the other hand, the 
conjugate variable, the charge, diminishes its intrinsic 
quantum fluctuations and can in principle be detected. 

The circuit we envisage is depicted in Fig. QJ. It con- 
sists of a three-junction loop frustrated by a magnetic 
flux close to half a flux quantum, following Ref. Q, one 
of the islands of which also serves as the island of a single- 
charge transistor. I.e., it is coupled to voltage sources by 
small Josephson junctions exhibiting Coulomb blockade, 
and is subject to an electrostatic gate voltage. For fix- 
ing the reference, the different parts of the device are 
connected to ground either capacitively or galvanically. 

In the limit we are interested in, both components are 
in their generic regimes, i.e. £j, qu bit > E c , qu bit E c , S et > 
-Ej.set- Here and henceforth, Ej and E c = 2e 2 /C are the 
Josephson and charging energies of the respective junc- 
tion. The interaction between flux and charge is via their 
mutual geometric phases. Charge moving around mag- 
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FIG. 1: Proposed circuit for the read out of a flux qubit (three 
junction loop on the left, junctions 1,2 and 3, controlled by 
the external flux $) at the degeneracy point using an SET 
(part on the right, other junctions). In order to establish 
the interaction, the qubit has to have (stray) capacitances to 
ground, C\, CV, and C$ as well. 



netic flux acquires an Aharonov-Bohm phase, a flux line 
moving around a charge acquires an Aharonov-Cashcr 
phas o 11 ! 12 . This analogy is a manifestation of vortex- 
charge duality. 

In order to accommodate loop constraints and Kirch- 
hoff rules, we start from a Lagrangian analysis. The 
Lagrangian of such a coupled circuit can be schemati- 
cally represented in terms of the junction phases (pi as 

£ {{4>i}, {</>?}) = £ch - £j - £s- The charging part 
2 -jr 

reads C = 4> where C is the capacitance 

matrix and the phase velocities have been arranged in 
a column vector. The Josephson part reads Cj = 
— J2j cos 4>i = Uj((f>i 1 (f>2,4 > 3) where the Josephson 
energies are given through the critical currents by Ejj = 
<f>o^c.j/27r. The work done by the voltage sources, finally, 

is counted as £5 = —^V T C g (j). In our case the capaci- 
tance matrix is, using the notation from Fig. 

C*8 + Ci + C3 — C\ Cg 
C = I — C\ C\ + C2 + Cj C-j 

C% Cj C7 + C$ + Cs 

(1) 

where Cs = 6*4 + C5 + Ce is the total capacitance of the 
transistor, and the source term is given by 






-C 5 V X - C 6 V T 



(2) 



From here, we can find the canonical variables and the 
Hamiltonian 
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The canonical 
now as H = H n 



momentum P is given by ~? = 
g . We want to rewrite this Hamiltonian 
I- -ff q ,sET + H S et- From Eq. we see 
that the coupling between the transistor and the qubit 
has the form of —id$ k P?,. This coupling term can, as 
well as the coupling to voltage sources, be transformed 
into phase factors by a gauge transformation ip(<j>) — > 
■0(0) exp i J Psdcjn . This results in a qubit Hamilto- 
nian in phase representation of the same shape as Ref. HH 
with a slightly modified capacitance matrix. Performing 

e A 



the two level approximation leads to H' = 



A* 



with A = |A|e lx where x = S<j>P^ an d S(f>, the separation 
of the minima in the potential landscape corresponding 
to the persistent current states^. In this gauge the charg- 
ing part of the transistor island has its standard form 
Ec(nL-7iR,-n g ) 2 + 2eV(nnL + (l-K)nR) , where n L / R 
is the number of charges that have tunneled through the 
left/right junction and rii — ur is the number of extra 
charges on the transistor island, k = Cr/(Cl + Cr) will 
be assumed to be 1/2 through the rest of this paper. 
C L = C 6 and Cr = C a + C 5 + C 6 + C 7 + C 8 are the ef- 
fective left/right capacitances of the transistor and n x is 
the dimensionless gate voltage. As -EcSET Ej.set, we 
can perform a two level approximation also for the tran- 
sistor island. In order to simplify the treatment of the 
external Josephson junctions, we observe that the phase 
X measures the charge tunneling through the SETs and 
obtain the qubit-transistor interaction: 



(H q + -f/q,SET)' - r 



Ae ~i<t>N 



Ae^ N 
— e 



(4) 



Here ^ is a constant set by the geometric capacitances 
of the circuit. At the flux degeneracy point, we have 
e = 0. By performing the unitary transformation U = 

exp (jj&z ® T z ^j , where a z is a Pauli matrix operating in 

in the qubit subspace and t z is a Pauli matrix operating 
in in the transistor subspace, we finally obtain 
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ch 



Ej S ET 1 ® COS( 



f+(Z + f) + f_(Z t + r?) 



iEjsET&x 55 sin < 



f+(l + f) -<?_(!* + ft) 
A- 



(5) 



Here, we already have assumed to operate at the flux 
degeneracy point e = and we have changed to the en- 
ergy basis such that a z now describes energy and a x flux 
cigenstatcs. n L / R count charges having tunneled to the 
left and right lead respectively and l/r transfer charges 
from the left/right lead to the island. Moreover, f± are 
the raising and lowering operators for the island charge. 
As we envision the readout to happen completely on the 
supercurrent branch (see below), we omit quasiparticle 
terms. For finding dissipative action and finally observ- 
ing a voltage signal, one has to include resistors described 
by the a and (3 operators, which are modeled as baths of 
harmonic oscillators along the lines of the familiar P(E) 
analysis of the Caldeira-Leggett modelii 

This is the starting point for our rate equation study 
of the readout process. The gauge transformation has 
turned Aharonov-Casher into Aharonov-Bohm phases. If 
the qubit is in one of the flux, a x , eigenstates, the envi- 
ronment will not induce any transitions of the qubit and 
the qubit state will only influence the phases of the two 
tunneling events. Such phase differences cannot be de- 
tected in the present geometry, as this would require a 
closed superconducting loop. Notably, qubit and island 
charge formally appear to be uncoupled, but are jointly 
affected by coupling to the superconducting leads. This 
is a result of our gauge transformation: In the previous 
gauge, the qubit eigenstates had different phases depend- 
ing on the eigenstate of f z , thus if the island charge de- 
scribed by t z flips the qubit can also make a transition. 
Note, that the qubit energy splitting has the same value 
in both t states. 

When Ej.set < E c .set and the dimensionless con- 
ductance of the junction in the transistor is small, we can 
treat the tunneling term in the Hamiltonian in Eq. (JjjJ 
perturbatively. The first two terms of Eq. JSJ) plus the 
free leads in the last two terms constitute the unper- 
turbed part. The device has to be biased at A = 5E C . 
The eigenlevels of the unperturbed system are shown in 
Fig. In the ideal case, the parameter cj) in Eq. (JSJ is 
exactly n/2. In this case tunneling through the junctions 
of the transistor is only possible if, at the same time, the 
qubit also switches. As a charge tunnels through first the 
left, then the right junction (corresponding to following 
the right-going arrow in Fig. then the left-going) a 
charge has been transfered through the transistor, while 
the qubit and the transistor returns to the original state. 
This lead to a quantum flip-flop process where the flips of 
the transistor island as current passes are accompanied 
by opposite transitions — flops — of the qubit: Charg- 
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FIG. 2: The eigenlevels of the unperturbed Hamiltonian. The 
arrows indicate the wanted transitions. The up/down arrows 
indicate the qubit states and the numbers 1/0 indicate the 
transistor states. 



ing the island relaxes the qubit, discharging the island 
excites the qubit. All other transitions are suppressed 
due to symmetry. Thus we get a current running through 
the transistor if the qubit is in the excited state, but no 
current in the ground state. By reading out the current 
running through the transistor the state of the qubit can 
be detected. As there is no mixing, we can take as much 
time to read out the state as we want. Remarkably, if af- 
ter the measurement the transistor island is uncharged, 
the post-measurement state is the qubit eigenstate corre- 
sponding to the measured observable without relaxation. 
This is the same outcome as in a QND measurement, 
even though the Hamiltonian does not satisfy the QND 
condition^ and the qubit does make transitions during 
the operation. This result is owed to the flip-flop sym- 
metry of the Hamiltonian. The tunneling probability per 
second is constant leading to a Poissonian current dis- 
tribution, which means that arbitrary large resolution is 
available in this ideal limit by reading out during arbi- 
trarily long time. 

An important note is that if the qubit is initially in 
a superposition of its eigenstates (remember, that eigen- 
states at the degeneracy point are themselves equal su- 
perpositions of current states), a single charge cycle in 
the transistor will destroy the coherence of the qubit. As 
a single charge cycle can in principle be detected with 
the method of Ref. Il5l this is a quantum-limited mea- 
surement. 

For the non-ideal case, the parameter (f> in eq. J3J) is not 
exactly n/2 and we have unwanted transitions that lead 
to mixing of the qubit. As </> depends on the capacitor 
setup, adjusting it with high precision appears to be dif- 
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ficult. In that case, there is a risk of making a transition 
to the | I 0) state by spontaneous emission to the en- 
vironment, see Fig. which halts the detection process 
and potentially falsifies the result. Thus, we need to find 
a detection scheme which suppresses spontaneous emis- 
sion, i.e. which has an effective spectral density peaked 
around zero. 

This can be done by using incoherent Cooper pair 
tunnelin g 1617 . By putting an impedance in series with 
the transistors the supercurrent of the transistor is sup- 
pressed. The fluctuations of the voltage across the junc- 
tion leads instead to an incoherent transfer of Cooper 
pairs. Still, the energy of the Cooper pair condensate de- 
fines a preferred tunneling window and leads to the de- 
sired band-pass property. The strong coupling this setup 
provides and the large qubit splitting compared to the 
temperature of the oscillator bath makes the bandpass 
properties quite accentuated which combined with the 
symmetry properties should result in a good signal to 
noise ratio. 

The incoherent rates can be calculated using standard 
P{E) theory for an Ohmic environment*^, using 

r E < -k B T/g 

P{E) = \ l/T \E\ < k B T/g (6) 
{2/gE E^k B T/g 

where g is the dimensionless conductance of the environ- 
ment. An example of P(E) is shown in Fig. [3J Thus, 
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FIG. 3: P(E) for an ohmic environment for different resistance 
of the environment. The junction capacitance used was 0.25 
fF and the temperature was 25 mK. Note the lin-log scale. 



the rates for the wanted processes are T ~ sin 2 (</>)/T 
whereas the unwanted one is r crror ~ cos 2 (0)2 / g E where 
E ~ 2E C is the energy gain of a transition to the dark 
|) state. This way, we can transfer a number of 
N = r/r orror = {gE/2T)tan 2 (4>) Cooper pairs through 
the device before an error occurs, setting an error win- 
dow for the measurement time of terror — N/ 1 where / is 
the average current. 

Using the results from Fig. OHwitli a resistance of 50 fi, 
this ratio can be estimated to several thousand Cooper 
pairs. Considering present day technology where single 
Cooper pairs can be counted^, this can be straightfor- 
wardly detected. 

This setup has several advantages compared to other 
setups. For instance it allows for a flux qubit read out on 
the flux degeneracy point as well a read out for arbitrary 
small qubit inductance. Because it uses geometric phases 
to couple the qubit and the pointer it does not require the 
qubit to move into a mixed charge-flux regime^. This 
should make the qubit insensitive to charge noise while 
at the same time maintaining insensitivity to flux noise 
to first order. The physics behind the coupling is also 
quite general meaning that it should be possible to adapt 
this read out scheme to any system containing conjugate 
variables. 

In this letter, we have demonstrated a novel and 
promising read out scheme for flux qubits using the con- 
jugate variable. It allows for read out on the flux degen- 
eracy point, leading to improved noise insensitivity. Us- 
ing incoherent Cooper pair tunneling as a form of band- 
pass filtering the two different qubit states can be dis- 
tinguished and detected. The strong coupling provided 
by the setup allows the bandpass filter to, in the low 
temperature regime, provide a good signal to noise ratio. 
We have identified an ideal working point in which the 
measurement has QND-like properties and is quantum- 
limited. 
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